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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

FLIGHT MEASUREMENTS OF THE DYNAMIC LATERAL 
AND LONGITUDINAL STABILITY OF THE BELL X-5 
RESEARCH AIRPLANE AT 58 . 7 ° SWEEPBACK 
By Edward N. Videan 

SUMMARY 


An investigation has been made of the dynamic stability of the 
Bell X-5 research airplane at 58*7° sweepback and at altitudes of 
UOjOOO feet and 25,000 feet over a Mach nvimber range of O. 5 O to 0.97* 

The results of this investigation show that the longitudinal oscil- 
latory motions are well damped over the entire Mach number range except 
for residual oscillations resulting principally from engine gyroscopic 
coupling with the lateral oscillatory mode. The lateral oscillatory mode 
is moderately well damped except for nonlinear damping characteristics 
above a Mach number of O. 8 O. The damping is high for large amplitudes 
but for sideslip angles of less than 2 ° the damping is poor. The air- 
plane is highly sensitive to small inadvertent control motions which 
produce apparently undamped small amplitude oscillations in the Mach 
number range above M = O. 8 O. 

By U. S. Military Specifications, the lateral oscillation is mar- 
ginal except in the nonlinear damping range where it is definitely 
unsatisfactory. 

The longitudinal and lateral frequency-response characteristics were 
obtained and some coupling effects were noted. Longitudinal frequency- 
response calculations were characterized by adjacently located double 
peaks in amplitude ratio. This behavior is satisfactorily explained by 
theoretical calculations of frequency response involving engine gyroscopic 
coupling to the lateral mode. Lateral frequency responses showed some 
dependence on the direction of the initial distiorbance. 
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INTRODUCTION 


The Bell X-5 research airplane was designed primarily for investiga- 
ting the effects of wing sweep angle on transonic airplane characteristic 
Accordingly, this airplane incorporates a wing which has sweepback vari- 
able in flight from 20° to 58 . 7 °. As part of the research program, 
fll^t measurements have been made with the X-5 airplane to determine 
the longitudinal and lateral dynamic stability characteristics over a 
Mach number range of O. 5 O to 0.9T sit an altitude of about 40,000 feet 
and with a wing sweep angle of 58 . 7 °. Flight measurements have also 
been made at 25,000 feet to determine altitude effects. This paper pre- 
sents the longitudinal and lateral dynamic characteristics in terms 
of P, T 1/2 altitudes of 40,000 feet and 25,000 feet. Some of the 
longitudinal stability derivatives are also presented for these altitudes 
The longitudinal and lateral frequency responses are given for the tests 
at an altitude of 40,000 feet. 


SYMBOLS 


b 


wing span, ft 


CPa = dCi/da 


lift coefficient, L/qS 


'e 


pitching-moment coefficient about airplane center 
of gravity 


- dCjji/da 
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“ C 

^rri0 " 

Cjn: rate of change of pitching moment due to engine gyroscopic 

^ moment with precessional yawing velocity 

normal-force coefficient 

CnA rate of change of yawing moment due to engine gyroscopic 

moment with precessional pitching velocity 

= dC^/d\|r 
*^rur, “ 

cycles for oscillation to damp to one-half amplitude 
cycles for oscillation to damp to one-tenth amplitude 
chord, ft 

mean aerodynamic chord, ft 

d/dt, differential operator with respect to time 

acceleration due to gravity, ft/sec^ 
pressure altitude, ft 

moment of inertia about Y-body axis through center of 
gravity, slug-ft^ 

moment of inertia about Z-body axis through center of 
gravity, slug-ft^ 

lift, lb 

Mach number 

mass, slugs 


Cl/2 

^l/lO 
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normal acceleration, g \mits 
period of oscillation, sec 

dynamic pressTjre, Ib/sq ft 

area of wing, sq ft 

time for oscillation to damp to one-half amplitude, sec 

time, sec 
velocity, ft/sec 

equivalent side velocity, 5^^ ft/sec 

angle of attack, deg or radians 
da/dt, radians/sec 
angle of sideslip, deg 

double amplitude of 3 in oscillatory mode 
single aileron position, deg 
elevator position, deg or radians 
rudder position, deg 
pitch angle, radians 

pitching angular velocity, de/dt, radians/sec 

air density, slugs/cu ft 

air density ratio 

phase angle, deg 

rolling angle, deg 

rolling angular velocity, radians/sec 
double amplitude of cp in oscillatory mode 
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angle of yaw, radians 

yawing angular velocity, d\|r/dt, radians/sec 
frequency, radians/sec 

IcpI ^7.3 

|p| VVa 


AIRPLANE AND INSTRUMENTATION 


A three-view drawing of the X-5 research airplane is shown in 
figure l(a), and a photograph is presented in figure l(b). Pertinent 
physical characteristics are described in table I. 

Standard NACA internal recording instruments were used to measure 
airspeed, altitude, normal acceleration, pitching, rolling, and yawing 
angular velocities, angle of attack, angle of sideslip, and control 
surface position. All recordings were synchronized by a common timer. 


The angular velocities were measiored by rate gyros with accuracies 
of 0.5 percent of instrument full scale. These instruments have the 
following ranges and dynamic characteristics: 


Angular velocity 


Range , 
radians /sec 


Natural frequency, 
cps (approx.) 


Damping 
(approx. ) 


Pitch ±0.5 
Yaw ±1.0 
Roll ±3.0 


9.0 

13.5 

17.0 


0.65 critical 
0.65 critical 
0.65 critical 


The flow direction recorders were vane-type pickups mounted on the nose 
boom. The undamped natural frequency is about 8 cycles per second and 
the damping ratio is approxi m ately 0 . 7 . 

Rudder, aileron, and elevator positions were meas\ired by trans- 
mitters linked directly to the control sinrfaces. The frequency response 
of the transmitter -recorder system has been measured and found to be 
flat to 20 cycles per second over the amplitude range of the control 
movements presented. 

Airspeed was measured from a calibrated nose-boom installation. 

The accuracy of the airspeed calibration is ±0.01 in Mach number. 
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TESTS 


The data were obtained from the transient oscillation of the air- 
plane in response to rudder pulses, elevator pulses, and aileron pulses. 
The rudder pulse magnitudes were 5° to 12^, the elevator pulse magni- 
tudes were 4° to 9^, and the aileron pulse magnitudes were 20^ to 28^. 
The time duration of all pulses was 0.5 second to 1.0 second. At all 
times the pilot attempted to have the pulse terminate with the control 
surface in its initial position and to hold the control surface fixed 
for the duration of the maneuver. Generally all other controls were 
also held constant throughout the maneuver. The pulses and resulting 
transient oscillations at pressure altitudes of 40,000 feet and at Mach 
numbers below 0.95 were obtained for initial 1 g trim flight conditions. 
To obtain data at speeds greater than M = 0.95 it was necessary to 
dive the airplane at angles up to 10^. However, these dive angles were 
s\iff iciently small so that the change in altitude would have negligible 
effect on the oscillation. A smaller number of pulses were made at 
25,000 feet to investigate altitude effects. Figure 2 presents the 
average normal-force coefficient as a function of Mach number for 

all pulses at both 25,000 feet and 40,000 feet. 


RESULTS AND DISCUSSION 
Longitudinal Stability 


Figure 3 presents representative time histories of the X-5 airplane 
transient -longitudinal response to an elevator pulse at 40,000 feet. 

These time histories were chosen to cover the important parts of the 
Mach number range of the flight test measurements. Figure 3 shows that 
the longitudinal motions immediately following the elevator pulse damp 
quickly over the entire Mach number range. However, the time histories 
also show that the longitudinal oscillations induce a lateral-directional 
mode of oscillation which in turn affects the initial transient- 
longitudinal oscillation. This coupling will be discussed later in 
conjunction with lateral stability and frequency response. 

Figure 4 presents the variations with Mach number of the period, 
time to damp to one-half amplitude, and cycles to damp to one-tenth 
amplitude for the longitudinal response of the airplane to an elevator 
pulse. These measurements were obtained from 0 over the first part 
of the transient oscillation immediately following the cessation of 
control movements. This portion of the curve was chosen for making the 
measurements in an attempt to minimize the effects of cross coupling 
with the lateral mode. It is impossible to exclude these effects 
entirely because, as the time histories show (fig. 5), the lateral- 
directional mode is disturbed almost as soon as the longitudinal mode 
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is distijrbed . This condition makes the measurement of T]_^2 particu- 
larly difficult and casts some doubt on the validity of the measurements 
of both period and time to damp to one-half amplitude, but it is consid- 
ered to be the best approximation to the true uncoupled longitudinal 
period and damping. The data at 40,000 feet (fig. 4) show that the 
period varies smoothly with Mach nxmiber, decreasing gradually from about 
2.0 seconds at M = 0.^6 to about 1.5 seconds at M = 0.92. Above this 
speed the period drops more sharply, reaching 1.1 seconds at M = 0.97* 

The data at 25,000 feet show that the period decreases by about 20 percent 
to 26 percent of the value at 40,000 feet. The approximate theoretical 
change with altitude is 29 percent, based upon the approximation that 
the ratio of the periods at the two altitudes is equal to the square root 
of the inverse ratio of dynamic pressure. 

The damping at 40,000 feet changes only slightly with Mach number, 
the time to damp to one-half amplitude varying from near 1.0 second 
at M = 0.55 to near 0.5 second at M = 0.97- This results in a vari- 
ation with Mach number of the number of cycles to damp to one-tenth 
amplitude as shown in figure 4. At 25,000 feet *^l/l0 

only small differences from the data at 40,000 feet. 

Figure 5 presents the variation of the longitudinal stability 
derivatives and with Mach nxmiber for the two test 

altitudes. These derivatives were calculated from the period and damping 
by the relations; 



The lift-c\irve slope used in the computation of Cjjj^ + was obtained 

from flight measurements as presented in reference 1. Reference 2 shows 
that above M = 0.6 and below the pitch-up boundary (ref. 5) the 
X-5 airplane possesses two static stability regions and a transitional 
region which are dependent upon Mach number and These regions 

(low lift - decreased stability, transitional stability, moderate lift - 
increased stability) are marked in figure 2. It may be seen that most of 
the pulse maneuvers performed at 40,000 feet are in or near the transi- 
tional region, whereas the maneuvers performed at 25,000 feet are in a 
region of constant though decreased stability. Values of Cu)^ (fig. '?) , 
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calculated from the data at 40,000 feet are, therefore, more subject to 
scatter than the data at 25,000 feet. Figure 5 shows that the agreement 
of Cjn^ calculations for the two altitudes is not perfect. It appeaxs 

that closest agreement is obtained between M = 0.8 and M = 0.9 where 
the average at the two altitudes is on the same side of the transi- 

tional region of stability (fig. 2). 

The values of Cjji^ + Cjq^, also presented in figirre 5, show some 

scatter which is common to this method of obtaining the derivative. The 
scatter is possibly aggravated by the aforementioned effects of coupling 
on the determination of T-Ly'2* However, the difference between the data 

at 40,000 feet and 25,000 feet is larger than the average scatter. 


Lateral Stability 

During the early flights at 58 . 7 ° sweepback, the pilots complained 
and records showed that the airplane oscillated almost continuously in 
the lateral-directional mode at small amplitudes, particularly at Mach 
nimibers above O. 8 O. An investigation was made to determine if this poor 
lateral dynamic behavior was affected by the magnitude of the disturbance 
and if the motions resulting from large disturbances were stable. 

Accordingly, rudder pulses of large amplitude were performed at 
40,000 feet over the Mach nimiber range from 0.52 to O .96 and rudder 
pulses of varying amplitude were performed at Mach numbers between 0.8 
and 0 . 9 . Throughout each pulse of the latter group the stick was held 
fixed by a mechanical stop to prevent inadvertent stick movements. 

Figure 6 presents typical time histories at 40,000 feet of the 
transient response of the X-5 airplane to rudder pulses over the Mach 
number range tested. In the measurement of damping from such time 
histories it was noted that ^\j2. large amplitude depended somewhat 

on the direction of the initial distiirbance. Also, though good damping 
as measured by usually present at large amplitudes, it was 

observed that the damping would decrease with decreasing amplitude 
resulting in near zero damping at small amplitudes (fig. 7)* To show 
the nonlinear effect, the damping was measured as shown by the hypo- 
thetical c\irve in figure 8. The "large amplitude" range and "small 
amplitude" range were separated at about ||3| =2°. The period and 
damping for these two portions of the oscillations are shown in figure 9. 
The characteristics of the residual oscillation are presented in fig- 
\rre 10. Figure 9 shows that the period varies smoothly with Mach number, 
ranging from 2.8 seconds at M = 0.52 to about 1.5 seconds at M = O. 96 . 
The damping measured over the large amplitude portion of the oscillation 
also shows a smooth and fairly consistent variation up to M « O. 75 . 
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For higher Mach numbers, despite scatter, there is a noticeable separation 
between values of T]_^2 large amplitude and small amplitude. This 

behavior appears to continue to the highest Mach numbers tested. The 
cycles to damp to one-half amplitude show a similar trend. In general, 
the damping of large amplitude lateral motion is good throughout the 
Mach nvimber range. 

The damping measured for small amplitudes (fig. 9) > however, shows 
appreciably large variations with Mach mrniber. The low amplitude damping 
begins to vary noticeably from the large amplitude damping at a Mach 
number of about 0.85 and reaches a maximum value about double the 

large amplitude T]_/2 M « 0.86 to M « 0.88. Above M » 0.95 the 

low amplitude damping increases rapidly and again generally approaches 
the large amplitude damping. The cycles required to damp generally 
follow a variation similar to the time required to damp. This variation 
between small amplitude and large amplitude damping thus defines a region 
of nonlinear damping extending from about M 0.80 to the test limit. 
However, the pilots report that the airplane apparently regains much of 
its damping in the small amplitude range above M = 0.93. 

Included as points in figure 9 are the results of the series of 
rudder pulses in the Mach number range of 0.80 to 0.90 in which rudder 
deflections of various amplitudes were made. As shown in the figure, 
points for large and small oscillatory amplitude as previously defined 
are fairly well grouped. The scatter existing in these two groups could 
not be Identified with the size of the rudder pulse. Therefore it is 
concluded that the magnitude of the rudder disturbance does not directly 
affect the damping to any extent, except as it affects the initial ampli- 
tude of the oscillations. 

Figure 10 presents the amplitude of the residual undamped oscil- 
lation of the lateral transient oscillation for the stick held fixed 
with the mechanical stop as well as with simple manual restraint. These 
data were obtained primarily at the end of the transient oscillations. 
However a few points, mostly those of the higher amplitude, are taken 
from undamped self -excited oscillations such as those of figure 7. These 
data show that with the stick restrained by the mechanical stop the 
residual oscillation amplitude reached a maximum of 0.5° as compared with 
a maximum residual amplitude in sideslip of about 2.5° when the stick is 
held manually. From these tests it would appear that small aileron move- 
ments are the major cause of the residual lateral oscillation. Figure 7 
does not appear to bear out this conclusion since all controls are fixed 
within the limits of the recorder as noted in figure 7 • However , as 
stated in the following pilots' comments, any small distiirbances such as 
gusts may initiate these oscillations. It will be noticed that the 
amplitude of the residual oscillation and the nonlinear damping character- 
istics follow the same variation with Mach number. Therefore, it would 
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appear that outside the nonlinear damping range, below M » 0.80, the 
higher damping at small amplitude of oscillation keeps the small disturb- 
ances from causing a continuous oscillation. 


Figure 11 presents the period and damping for the response to rudder 
pulses at 25,000 feet. It may be noted from the figure that the cycles 
to damp to one -half amplitude show about the same variations as shown by 
the large amplitude points at 4-0,000 feet. However, at the lower alti- 
tudes and lower (fig- 2) the nonlinear damping effects have disap- 

peared and the damping remains constant with amplitude over the whole 
amplitude range. 

The data shown in figures 9 to 11 are in substantial agreement with 
the pilots' comments on the dynamic stability of the X-5 airplane. NACA 
pilots who have flown the X-5 airplane feel that the Mach number range 
of 0.85 to 0.92 is dynamically less stable than at speeds above this 
range. The pilots agree that in the least stable region the airplane is 
subject to almost continuous lateral oscillations of low amplitude with 
zero damping, similar to the time history of figure 7 • The pilots report 
that gusts, changes in power setting, or small control motions initiate 
these undamped oscillations. 


Figures 12 and 13 present the new Military Specifications for 
dynamic lateral stability (ref. 4) and the recently superseded U. S. Air 
Force criteria for dynamic lateral stability (ref. 5), respectively. The 
new Specifications (fig. 12) relate the reciprocal of the nimiber of cycles 
to damp to one-half amplitude to the ratio of roll angle to side velocity. 


which is given by 


ItI 


M 

131 


57.5 


The value of cp was obtained by 


|ve| I PI VVa 

Integrating the roll velocity, whereas p was measured directly from a 
flow direction recorder. The points presented are calculated for the 
various damping regions of the transient oscillation previously mentioned. 
According to the new Specifications, some points fall into the satis- 
factory region and others into the unsatisfactory region. At the lowest 
Mach number range the points fall in the satisfactory region but are 
near the border of the ^satisfactory region. In the other speed ranges 
the majority of points fall into the unsatisfactory region. 


The superseded Air Force criteria (ref. 5) is presented in figure 13 . 
A representative number of points over the Mach nmber range have been 
plotted in this figure. The majority of points fall into the vinsatis- 
factory region but are fairly close to the boundary between the unsatis- 
factory and satisfactory region. 


It will be noticed that by both the old and new criteria the dynamic 
lateral stability of the X-5 airplane is predominantly unsatisfactory. 
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In evaluating the airplane, the pilots feel that the lateral dynamic 
stability characteristics are unsatisfactory, but tolerable, except in 
that portion of the speed range where the damping becomes noticeably 
nonlinear . In this region the low damping at small amplitude coupled 
with a large roll-to-sideslip ratio makes the behavior intolerable. 


Coupling of Oscillatory Motions 

As shown by the time-history plots (figs. 3^ 6, and 7)^ there is 
coupling between the lateral and longitudinal oscillatory modes. In 
general, coupling can be caused by inertial coupling effects, aerodynamic 
coupling, or engine gyroscopic effects. Although the relative magnitude 
of each type of coupling has not been determined for the X-5 airplane, 
examination of the time histories would indicate that both gyroscopic 
and aerodynamic coupling are present in appreciable amounts. 

Figures 5(b) and 5(c) show the coupling from longitudinal to lateral 
directional mode for an elevator pulse. Note that an up-elevator pulse 
produces a right yawing motion, whereas a down-elevator pulse produces 
a left yaw. Since the aerodynamic coupling to the lateral mode from a 
disturbance in pitch is usually negligible, this coupling is believed 
due to gyroscopic effects. The initial directions of coupled motion for 
the two disturbance directions agree with those expected from the gyro- 
scopic processional torque of the jet engine. 

Figures 6(a) and 6(b) show the coupling from the lateral mode to 
the longitudinal mode for rudder pulses. These time histories show that 
a right rudder pulse produces a small initial decrease in angle of attack 
and a left rudder pulse produces a large increase in angle of attack. 
Again gyroscopic coupling is implied; however, there are large initial 
amplitude and phase differences in the coupled oscillation for right 
and left rudder pulses, suggesting the probability of both aerodynamic 
and gyroscopic coupling. Sideslips performed with the airplane have 
shown no appreciable pitching moments caused by sideslip, indicating no 
static aerodynamic coupling. 

In the X-5 airplane configuration the mass is concentrated to a 
large extent in the fuselage and near the center of gravity. This phys- 
ical characteristic makes the moments of inertia about the Y- and Z-axes 
relatively small as compared with a conventional jet-powered airplane of 
comparable weight (determined experimentally by the method of ref. 6, 

ly = slug-ft^, I2, = 10,110 slug-ft^). Therefore, engine gyro- 

scopic torques are effective in producing disturbances in pitch and yaw. 
It will also be noted by comparing figure 4 and figure 9 that the periods 
of the longitudinal and lateral oscillations are almost equal. This con- 
dition will produce a near resonant state for the gyroscopic coupling 
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torque and resultant coupled oscillation. In this case a poorly damped 
oscillation in either the longitudinal or lateral mode voiild produce a 
forced resonant oscillation in the other mode. For example, an elevator 
pulse in the directional nonlinear damping region would couple a disturb- 
ance to the lateral mode, producing a small poorly damped oscillation 
which is coupled back to the longitudinal mode and appears there as a 
near resonant forced oscillation. For these reasons engine gyroscopic 
effects are considered important in an analysis of the characteristics 
of the X-5 airplane. 


Frequency Response 

The transient responses of the X-5 airplane have been analyzed by 
means of the Fourier transform and the frequency response obtained. 

These calculations were made in order to present the dynamic character- 
istics of the airplane in a form that would be more usable to automatic 
control system designers. The calculations were made on an IBM machine 
calculator by a numerical integration procedure using the methods of 
reference 7 . The data selected for analysis were from the same tran- 
sient responses previously discussed and were analyzed about the airplane 
body axis. In addition to the above transients, a. series of aileron 
pulses were made over the Mach mmiber range for frequency response only. 
The airplane responses analyzed and presented in this paper are &/5e^ 

^/^a- 

Reference 8 shows that the frequency-response characteristics are 
invariant with control pulse shape and size only when the system is 
linear. It has been shown previously that the lateral mode of the 
X-5 airplane is nonlinear in damping over part of the amplitude range 
of the oscillation. The frequency-response curves therefore represent 
the airplane only for the range of rudder amplitudes used in this 
investigation. 

Presented in figure l4(a) are typical results of the investigation 
in the longitudinal mode, showing the frequency response of the pitching 
velocity for an elevator input S/bg as obtained by analyzing an elevator 
pulse maneuver. This figure shows the amount of scatter obtained in the 
calcTulated response points and the fairing made in such cases. Fig- 
ures l4(b) to lU(d) present composites of these curves showing the Mach 
number variation. The double peak found in these responses is rather 
unusual and at first seems to indicate the presence of two modes of 
oscillation with natirral frequencies fairly close together. The simi- 
larity in magnitude of the natural frequencies of the longitudinal and 
directional modes plus the large mode coupling appears to offer a basis 
for explanation of the shape of the longitudinal frequency-response 
cvirves. However, the time histories of pitching velocity do not appear 
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to contain two frequency modes, nor is there any natural frequency in 
pitch or yaw closely corresponding to the lower frequency-response peak. 
In an effort to show the effects of engine gyroscopic coupling, tests 
were made performing elevator pulses at a given speed and altitude with 
the engine idling and at maximum revolutions per minute. To achieve as 
large a difference in engine speed as possible, the tests were performed 
at 15,000 feet where the engine can be idled as low as 60 percent of 
maximum rpm. The tests consisted of a pulse made with the engine idling 
at 60 percent of maximum rpm while the airplane speed was decreasing 
through M = 0.6. Another pulse at 100 percent rpm was made while the 
airplane speed was increasing through M = 0.6. The frequency-response 
calculations for these pulses are presented in figure 15 . The frequency- 
response plot shows that the difference between the maximum and minimum 
in the double peak is slightly less for the low rpm than for the 100 per- 
cent case. This difference wo-uld tend to substantiate the fact that 
there is an engine gyroscopic moment acting on the airplane. However, 
considering the accuracy of the tests, the differences are too small for 
any definite conclusions to be drawn. 

In a further attempt to determine the cause of the double peak 
longitudinal frequency response, calculations of longitudinal frequency 
response were made from theoretical transfer functions derived from 
simplified equations of motion as shown in the appendix. Two sets of 
calculations were made, assuming first a two-degree-of -freedom system 
and then a three-degree-of -freedom system. The two -degree -of -freedom 
system consisted of the classical longitudinal equations and assiomed 
constant forward velocity. The three-degree-of -freedom system included 
the same two previous equations, modified to include a moment term due 
to gyroscopic coupling from yawing velocity. The third degree of freedom 
was assumed to be yaw about the Z-axis. This equation also included a 
moment term due to gyroscopic coupling from pitching velocity. The 
transfer function 0/be was obtained for both the two- and three-degree- 
of -freedom cases and the frequency response was calculated by using 
longitudinal and simplified lateral derivatives which were obtained by 
assuming only a single degree of lateral freedom. These lateral deriv- 
atives were calculated from the flight data as described in the appendix. 
These theoretical frequency-response calculations were made for two 
Mach numbers, M = O.7O and M = O.87, at 40,000 feet and are presented 
in figures l6(a) and l6 (b). The responses are admittedly rough approxi- 
mations to the actual airplane frequency responses, since at least 
two degrees of lateral freedom are neglected along with all inertial 
and aerodynamic coupling terms. However, it is evident that the dominant 
airplane response characteristics are demonstrated by the simplified 
three-degree-of -freedom calculations. Figure l6(b) shows that the fre- 
quency responses obtained from the transfer function of the three-degree- 
of -freedom system match in general shape and amplitude the frequency 
responses obtained experimentally (fig. l4). The agreement of frequency 
location for the corresponding peaks in these figures is not too good. 
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Results of calculations of ^&e "the classical two-degree-of-freedom 

system, neglecting coupling, are given in figure l 6 (a). It is interesting 
to note, in comparing the two-degree-of-freedom system to the three- 
degree-of -freedom system, that the frequency locations of all the peaks 
differ. The transfer functions were factored to determine the roots of 
the numerator and denominator. The classical two-degree-of-freedom system 
was found to be of the form 

A. = (D + 7) 

(d2 + ad + b) 

where the denominator has a pair of conjugate complex roots. The three- 
degree-of -freedom system transfer function was found to have the form 

■ (D . ri) (b" ^ a^d ^ 

^d2 4- A 2 D + + A 3 D + 

with a real root and a pair of conjugate complex roots in the nvunerator 
and two pairs of conjugate complex roots in the denominator. The complex 
roots indicate that this system has two natural frequencies. However, 
as shown by the values of the roots given in the appendix, the conjugate 
complex pair in the numerator is close to the value of one of the conju- 
gate complex pairs in the denominator and tends to cancel the effects of 
the latter. 

Figures I 7 , I 8 , 19, and 20 present, respectively, the results of 
the experimental frequency-response calculations for rudder pulse 

ijfyfej., and Since some differences were noted on the time 

histories between the airplane response to right and left piilses , the 
frequency-response characteristics are presented separately for right 
and left pulses. Also presented in each figure is a plot showing the 
individual calculated points for a specific type of response. Plots 
of p/5r presented in figures 17 (a) to 17 (c). It will be noticed 

that the frequency-response variations with Mach number are not always 
consistent. It is not known whether this lack of consistency is due to 
nonlinearity in some aerodynamic parameters or to scatter in the data. 

The composites for right and left rudder pulses show little difference 
except that the left pulses have a smaller amplitude ratio peak at the 
lowest Mach numbers . The general remarks made about P/^r true 

also for and (figs. I 8 and 19). However, for the rolling 

velocity response (fig. 19 ) the right pulse amplitude ratio peaks 

appear to increase somewhat over the left at the highest Mach numbers. 

The pitching response to rudder disturbance presented in 
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figure 20. Although the differences in amplitude ratio are small for 
right and left pulses, the differences in phase angle are large. These 
differences were noted previously from an examination of time histories. 

Aileron pulses in both directions were also made and responses were 
analyzed in frequency response fom. The results for rolling velocity 
response are presented in figure 21. Again, large differences 

are noted in the curves for right and left pulse groups. An effort was 
made to obtain the frequency response results were not of 

sufficient quality to present. 

Because of the complicated nature of the transfer functions of this 
airplane and the presence of some nonlinear derivatives, no attempt has 
been made to determine the lateral stability derivatives, except some 
simplified derivatives for the longitudinal frequency-response calcu- 
lation described in the appendix. 


CONCLUDING REMARKS 


An investigation of the dynamic stability of the X-5 research air- 
plane at 58 * 7 ^ sweepback at altitudes of 40,000 feet and 25,000 feet over 
a Mach number range of 0.50 to 0.97 shows the following: The longitu- 

dinal motions are well damped over the entire Mach mmiber range tested 
except for residual oscillations resulting principally from engine gyro- 
scopic coupling with the lateral oscillatory mode. This engine gyroscopic 
coupling results in motion in both the longitudinal and lateral oscilla- 
tory modes for either a longitudinal or lateral disturbance. The lateral 
oscillatory mode exhibits moderately good damping except for nonlinear 
damping characteristics above a Mach number of 0.80 where the motion is 
well damped at large amplitudes and poorly damped for double amplitude 
sideslip angles smaller than 2^. Small inadvertent aileron control 
motions often produce apparently undamped small amplitude oscillations 
in the Mach number range above 0.80. 

^ the current military aircraft dynamic stability requirements, 
the lateral oscillation is unsatisfactory over most of the Mach number 
range. The pilots concinr with this finding and report the airplane to 
be particularly intolerable in the range with nonlinear damping 
characteristics . 

The longitudinal frequency response S/bg has a large and unusual 
double peak near the short period natural frequency. Theoretical calcu- 
lations have shown that engine gyroscopic coupling effects are princi- 
pally responsible for the shape of the frequency-response curves. 
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Lateral frequency-response calculations show some differences for right 
and left rudder or aileron pulse disturbances, particiilarly the longi- 
tudinal response to a rudder pulse 0/6r and the lateral response to 
an aileron pulse tp/Sa* 


High-Speed Flight Station, 

National Advisory Committee for Aeronautics, 
Edwards, Calif., Avigust 5, 1955* 
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APPENDIX 


DERIVATION OF A SIMPLIFIED THEORETICAL TRANSFER FUNCTION 
DESCRIBING THE ENGINE GYROSCOPIC COUPLING EFFECTS 


The form of the equations of the tvo-degree-of -freedom system and 
its trsinsfer function will not be presented in this paper since it is 
generally well known and may be found in reference 6. All the derivatives 
and airplane constants used in the calculations are contained in the 
derivatives and constants listed for the three-degree-of -freedom system 
(table II) . 

The equations of motion assumed for the three-degree-of -freedom 
system are: 


f B(e - o.) . Ci^a (1) 

D^0 = + Cjn^Da + CnigD© + (2) 

^ = Cn^^lr + Cn^D^r + Cn§D0 (5) 

These equations assume that the stability axes coincide with the airplane 
principal axes of inertia and that there are two-degrees-of -longitudinal 
freedom and one-degree-of-lateral freedom. In formulating these equa- 
tions it was assumed that the rotating engine mass produced gyroscopic 
moment proportional to the precessional angular velocity where 

p 

R = axial moment of inertia, slug-ft"^ 

G = spin of engine mass, radians/sec 
A = precessional angular velocity, radians/sec 
Q = moment, ft -lb 

The coupling derivatives and were obtained as follows: Let 

M = pitching moment 
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then 



letting 


M = Q = RGA 




= mL 

qSC 


dC, 


m RG 
— 


d(\^) qSc 


where 


A = ijr 


The yawing-moment derivative Cng is obtained in a similar manner so 
RG 

that CjiA = — The derivative corresponds to the usual deriv- 


"“e qSb • 




ative C 


np- 


In assimiing a single-degree-of-lateral freedom -p = i|f 




For this simplified case 


-I. 


qSb 


^ 'o.69? \g 



This relationship was used to calculate Cn^. The other lateral deriv- 
ative Cn^ was found to have a very small effect and was estimated. 

The other derivatives in table II were obtained from figures and refer- 
ences presented in this paper. In solving the equation, let 


qSb ^ 
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Divide both sides of eqiiation (l) by equation (2) by k£, and equa- 

tion (3) by kj. By solving these equations S/fee Is obtained as 


^ Aj_D^ + A2D^ + AjD + A4 

D^ + A^D^ + AgD^ + A^D + Ag 


^here 


At = 


‘"nu 

^ ^e 


1 k^ ^1^2 


Ag - 


■% 



Cl 

1 

o 

kqkg 

kj^kgk^ 

kik2 

k2kj 




A5 = - 








,kgk3 ^ 


AU - - [^ct^n^CLg^ + Cn^CL^Cj^ J 


Ar = 


^0 


kc 


Ct C 

^5_ 


Cl Cm, c^. 

a 6 . ^ 


CnrCL '^n. 


CmrCnA Cm 


_ g, o - V ¥ g _ 

^ ~ k-^k^ kgkj - k^kgkj ■ k 


L I xj AJ. I ■‘•AA r ii A — 

\|r g ^ \lr 6 a 


^2^3 


A>7 - 




kok- 


2^5 




k2k3 


C^ Ct C^. C^ C^ 

^8 " k^k^k^ kg^3 
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The frequency response was obtained from this transfer fimction by substi- 
tuting into the transfer function the proper values of the constants and 
D = icu. 

Substituting the aerodynamic coefficients at M = O.87 into the 
transfer f\anctions and factoring yields the following expressions; Two 
degrees of freedom 


-1^.082D + 1.^10 


D^ + 2.852OD + 19.3572 


Three degrees of freedom 


15.082D^ - 8.8170^ - I95.O82D - 19.451 


D^ + 3.557D^ + 55.5670^ + 47.009D + 249.1+46 


or 



(d^ + 2.56D + 21.8 ) (d^ + 0.80D + 11.3) 
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TABLE I 

PHYSICAL CHARACTERISTICS OF BELL X-5 AIRPLAME AT A SWEEP ANGLE OF 53.7° 


Airplane: 

Weight, lb: 

Full fuel 10,006 

Less fuel 7)89^ 

Power plant: 

Axial-flow turbojet engine J35-A-17 

Center-of -gravity position, percent M.A.C.: 

Full fuel ^5.0 

Less fuel ^5-5 

Moments of inertia for 58 . 7 ^ sweep (clean configuration full fuel, 

p 

body ajcis), slug-ft : 

About X-axis 5^185 

About Y-axis 

About Z-axis 10,110 

Overall height, ft 12.2 

Overall length, ft 35*6 


Wing: 

Airfoil section (perpendicular to 38.02-percent-chord line): 

Root 

Tip 

Sweep angle at 0.25 chord, deg - 

Area, sq ft 

Span, ft 

Span between equivalent tips, ft 

Aspect ratio 

Taper ratio 

Mean aerodynamic chord, ft 

Location of leading edge of mean aerodynamic chord, fuselage station 

Incidence root chord, deg 

Dihedral, deg 

Geometric twist, deg 

Wing flaps (split): 

Area, sq ft 

Span, parallel to hinge center line, ft 

Chord, parallel to line of symmetry at 20° sweepback, in.: 

Root 

Tip 

Travel, deg 

Slats (leading edge divided): 

Area, sq ft 

Span, parallel to leading edge, ft 

Chord, perpendicular to leading edge, in.: 

Root 

Tip 

Travel, percent wing chord: 

Forward 

Down 

Aileron (^5 percent internal-seal pressure balance): 

Area (each aileron behind hinge line), sq ft 

Span parallel to hinge center line, ft 

Travel, deg 

Chord, percent wing chord 

Moment area rearward of hinge line (total) , in.5 


, . NACA 614-^^Q^AOll 
NACA 61 +^qq^A008.28 

58.7 

183.7 

20.1 

19.3 

2.20 

0.411 

9.95 

101.2 

0 

0 

0 

15.9 

6.53 

30.8 

19.2 

60 

14.6 

10.3 

11.1 

6.6 


10 

5 

3.62 

3.15 

tl5 

19.7 

4,380 
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TABLE I.- Concluded 

PHYSICAL CHARACTERISTICS OF BELL X-5 AIRPLANE AT A SWEEP ANGLE OF 58 . 7 ° 
Horizontal tail: 

Airfoil section (parallel to fuselage center line) . . . NACA 65 AOO 6 

Area, sq ft 31-5 

Span, ft 9*56 

Aspect ratio 2.9 

Taper ratio 0.371 

Sweep angle at 0.25-percent chord, deg 45 

Mean aerodynamic chord, in 42.8 

Position of 0.25 mean aerodynamic chord, fuselage station . . 355*6 

Stabilizer travel, (power actuated), deg: 

Leading edge up 4.5 

Leading edge down 7*5 

Elevator (20.8 percent overhang balance, 31.5 percent span); 

Area rearward of hinge line, sq ft 6.9 

Travel from stabilizer, deg: 

Up 25 

Down 20 

Chord, percent horizontal tail chord 30 

Moment area rearward of hinge line (total), in. 5 4,200 

Vertical tail: 

Airfoil section (parallel to rear fuselage center 

line) NACA 65 AOO 6 

Area, (above rear fuselage center line), sq ft 25.8 

Span, perpendicular to rear fuselage center line, ft 6.17 

Aspect ratio 1.47 

Sweep angle of leading edge, deg 46.6 

Fin: 

Area, sq ft 24.8 

Rudder (23.1 percent overhang balance, 26.3 percent spein): 

Area rearward of hinge line, sq ft 4.7 

Span, ft 4.43 

Travel, deg t35 

Chord, percent horizontal tail chord 22.7 

Moment area rearward of hinge line, in. 5 3,585 
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TABLE II 

PARAMETERS USED IN ESTIMATING AIRPLANE FREQUENCY RESPONSE 

[hp = 40,000 feetj 


M 

0.70 

0.87 

V, fps 

680 

845 

q, Ib/sq ft 

135 

199 

S, sq ft 

185.7 

185.7 

c, ft 

9.95 

9.95 

b, ft 

20.1 

20.1 

radian"^ 

2.80 

2.865 

radian”^ 

-1.976 

- 1.848 

radian”^ 

-.4870 

-.4854 

Cm^, (radieins/sec)"^ 

-.0070 

-.0070 

Cuig, radian"^ 

-.3438 

-.3780 

(radians/sec)"^ 

-.0592 

-.0597 

CmA^ (radians/sec )“^ 

-.0522 

-.0527 

Cjn^, (radians/sec)"^ 

-.0556 

-.0372 

Cng, (radian/sec ) “^ 

.0278 

-.0187 

radian -1 

-.2550 

-.1862 

Cn^, (radians/sec)-l 

-.0100 

-.007 
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(a) Three -view drawing. 

Figure 1.- Bell X-5 research airplane at 56*7° sweepback. 
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(b) Photograph. 
Figure 1.- Concluded. 


L-87906 
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Figure 2.- Average Cjj during pulse maneuver. 


CONFIDENTIAL 



28 


CONFIDENTIAL 


NACA RM H55H10 



Time, t, sec 


(a) M = 0.70. 

Figxxre 5*~ Time history of airplane transient oscillation in response to 

elevator pulse at 40,000 feet. 
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(b) M = 0.85. 

Figure 5-- Continued. 
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Time, t, sec 

(c) M = 0.96. 
Figure Concluded. 
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Figure 4.- Period and damping for the longitudinal response of the 
X-5 airplane to an elevator pulse. 
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Figure 5*- Variation of the longitudinal stability derivatives 

and + C m^ with Mach number. 
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(a) M = 0 . 70 . 

Figure 6.- Time histories of transient oscillations in response to a 

rudder pulse at 40,000 feet. 
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(b) M = 0.67. 
Figvire 6.- Continued. 
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(c) M = 0 . 89 . 
Figure 6 .- Continued. 
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(d) M = 0 . 89 . 
Figure 6 .- Continued. 
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(e) M = 0.9ii-. 
Figure 6.- Continued. 
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(f) M = 0 . 96 . 
Flgiore 6 .- Continued. 
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(g) M = 0 . 96 . 
Figure 6.- Concluded. 
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Figure 7*“ Small amplitude undaji 5 )ed oscillation with all controls held 
fixed to less than 0.05^ motion. M = 0.88; hp = ^0,000 feet. 
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Amplitude of oscillation 





Figure 8.- Transient oscillation indicating various damping regions. 
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( a ) Large amplitude . 


(b) Small amplitude. 


Figure 9.- Period and damping for the lateral response of the X-5 airplane to a rudder pulse input 

at 40,000 feet. 
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Figure 10.- Variation with Mach number of the amplitude of the residual 

lateral oscillation at ^0^000 feet. 
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Figure 11.- Period and damping for the lateral response of the X-5 air- 
plane to a rudder pulse input at 25^000 feet. 
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Figure 12.- Comparison of the dynamic lateral stability of the X-5 air- 
plane with Military Specifications, reference k. 
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FigLire 15 •- Comparison of the lateral dynamic stability of the X-5 air- 
plane with superseded Air Force criteria, reference 5. 
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(a) Sample calculation showing calculated frequency response points. 

M = 0.69; hp = 40,000 feet. 

Figure l4.- Longitudinal frequency response as calculated from 
longitudinal transient responses. 
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(L) 


Figure l4.- Continued. 
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(c) 


Figure Continued. 
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(d) 

Figure l4.- Concluded. 
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Figmre 15*- Longitudinal frequency response with engine idling and at 
maximum rpm at 15,000 feet and Mach nimber of 0.6. 
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(a) Two degree-of -freedom system. 

Figiire l6.- Longitudinal frequency response calculated from theoretical 
transfer function, hp = 40,000 feet. 
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(b) Three-degree-of-freedom system. 
Figure l6.- Concluded. 
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(a) Sample calculation showing calculated frequency response points 
at M = 0.75; hp = 40,000 feet. 

Figure 17*- Lateral frequency response P/Sj, as calculated from 

lateral transient responses. 
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(b) Frequency response for right rudder pulses. 
Figure 17«- Continued. 
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(c) Frequency response for left rudder pulses. 
Figure 17 •- Concluded. 
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(a) Sample calculat.ion showing calc\ilated frequency response points 
at M = 0.75; hp = 40,000 feet. 

Figure l8.- Lateral frequency response as calculated from 

lateral transient responses. 
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(b) Frequency response of right rudder pulses. 
Figure l8.- Continued. 
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(c) Frequency response of left rudder pulses. 
Figure l8.- Concluded. 
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(a) Sample calculation showing calculated frequency response points 
at M = 0.75; hp = 40,000 feet. 


Figure I 9 .- Lateral frequency response 9 / 63 , as calculated from 

lateral transient responses. 


CONFIDENTIAL 


NACA RM H55H10 


CONFIDENTIAL 


6l 




(.d; Frequency response for right rudder pulses. 
Figure I 9 .- Continued. 
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(c) Frequency response for left rudder pulses. 
Figure I9,- Concluded. 
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(a) Sample calculation showing calculated frequency response points 
at M = 0.75; hp = ^0,000 feet. 

Figure 20.- Longitudinal frequency response as calculated 

from lateral transient responses. 
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(b) Frequency response for right rudder pulses. 
Figure 20.- Continued. 
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(c) Frequency response for left rudder pulses. 
Figure 20.- Concluded. 
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(a) Sample calculation showing calculated frequency response points 

at M = 0.88; hp = 40,000 feet. 

Figure 21.- Lateral frequency response cp/6g^ as calculated from 

lateral transient responses. 
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(b) Frequency response for right aileron pulses. 
Figure 21.- Continued. 
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(c) Frequency response for left aileron pulses. 
Figure 21.- Concluded. 
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